Chiral Mixed Phase in Disordered 3d Heisenberg Models 
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Using Monte Carlo simulations, we compute the spin stiffness of a site-random 3d Heisenberg 
model with competing ferromagnetic and antiferromagnetic interactions. Our results for the pure 
limit yield values of the the critical temperature and the critical exponent v in excellent agreement 
with previous high precision studies. In the disordered case, a mixed "chiral" phase is found which 
may be in the same universality class as 3d Heisenberg spin glasses. 
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The critical behavior of magnetic systems with 
quenched disorder has been of considerable theoretical 
and experimental interest for over 30 years. Disorder can 
lead to a competition between two order parameters and 
to a phase diagram which has regions in which each one 
orders independently as well as a mixed phase where both 
order simultaneously. These phases often meet at multi- 
critical points and can have various properties which de- 
pend on the spin dimensionality. Aharonji has recently 
reviewed some of the old and new results on multicritical 
points with special attention to high-T c materials. How- 
ever, the role played by quenched randomness remains a 
open question. 

Early renormalization group studies of systems with 
isotropic interactions and quenched disorder indicated 
that competing phases existed with multicritical points 
described by complex exponents^ or the absence of a sta- 
ble fixed poinl£. Various scenarios were suggested which 
included first order transitions, smeared transitions or 
spin glass ordering. For many years it was believed that 
a spin glass phase does not exist in 3d XY and Heisen- 
berg systems. However, numerical studies* during the 
last decade have suggested that chiral degrees of free- 
dom are important in spin glasses which have contin- 
uous degrees of freedom. More recent work££ indicates 
that Heisenberg systems exhibit a finite temperature spin 
glass transition in three dimensions where both the spin 
and chiral degrees of freedom order simultaneously. The 
correlation length critical exponent is estimated to be 
v = 1.1(2) which differs substantially from that of the 
pure 3d Heisenberg ferromagnet. 

In the present work we report a finite size scaling study 
of the spin stiffness of a 3d site random isotropic Heisen- 
berg model introduced previously by Matsubara et. al7 . 
The model describes a mixture of A and B magnetic ions 
randomly distributed at the sites of a simple cubic lattice 
with concentrations p and 1 — p respectively. The ex- 
change bonds between neighboring ion pairs are defined 
so that + J, —J, and +J correspond to A — A, B — B, 
and A — B (or B — A) pairs respectively. The Hamitonian 
can be written as 



where Si = (Sf, Sf, Sf) represents a classical three com- 
ponent spin of unit magnitude located at each site i and 
£,: = 1 or — 1 for A or B ions respectively. The average 
value of £j over the lattice is < e% >= 2p— 1. As discussed 
by the previous authors, the Hamiltonian has a symmetry 
with respect to p and (1 — p) and hence the phase dia- 
gram is symmetric about p = 0.5 as shown schematically 
in Fig. 1 of their paper. There are three ordered phases 
corresponding to conventional ferromagnetism, antifer- 
romagnetism and a mixed intermediate phase in which 
both a longitudinal ferromagnetism and transverse anti- 
ferromagnetism exist. The previous work used a finite 
size scaling analysis of the order parameter cumulants to 
determine the critical behavior at these transitions. 

We show that the spin stiffness for vector spin models 
can provide an independent alternative to the usual or- 
der parameter cumulants for locating the critical temper- 
ature and determining the critical exponent v associated 
with the correlation length. Our results for the phase di- 
agram and critical exponent v differ from those reported 
previously. 

We divide the bipartite lattice into two identical inter- 
penetrating sublattices a and (3 and consider the follow- 
ing order parameters 
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where rfi is the total magnetization per site and rh s is 
the staggered magnetization per site. These two order 
parameters describe the pure non-disordered phases at 
p = 1, respectively. Due to the symmetry about p = 0.5 
we shall only consider the region p > 0.5. 

The spin stiffness (helicity) tensor can be written 8 
as an equilibrium correlation function by taking second 
derivatives of the free energy with respect to the strength 
of an imposed spin gradient involving a twist about a par- 
ticular direction in spin space. We choose an orthogonal 
set of spin axes which correspond to unit vectors in the 
directions of m s x to, rh x (to s x to) and to which we de- 
note as 1, 2 and 3 respectively. Hence 2 is a unit vector in 
the direction of the component of the staggered magne- 



tization which is perpendicular to m. The spin stiffness elements can be written as 
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and the angular brackets 



< ... > indicate a thermal average. The superscripts 
a, b denote the components of the spin vectors in the 
orthonormal system 1, 2, 3. The ey are unit vectors along 
neighboring bonds and u is the spatial direction of the 
gradient in the lattice. 

Pure Heisenberg Model 

We first consider the pure limit p = 1 which describes 
the 3d ferromagnetic Heisenberg model. In this case 
there is no staggered magnetization and the two direc- 
tions perpendicular to m are equivalent. We find that 
the spin stiffness tensor is diagonal with the principal 
values p\ = P2 7^ /°3- The three principal stiffnesses are 
computed for simple cubic lattices with periodic bound- 
ary conditions and linear dimensions L ranging from 6 
to 32 using a single-spin flip heat bath algorithm 9 . We 
used 1 — 5 x 10 5 Monte Carlo steps (MCS) for performing 
our measurements after discarding the first 5 x 10 4 MCS 
to reach thermal equilibrium. In order to have continu- 
ous functions of temperature we also applied histogram 
reweighting techniques 1 ^ near T c using longer runs of 
(2 — 5 x 10 7 ) MCS in order to improve the statistics. Since 
we have continuous degrees of freedom, the histograms 



4.5 r 

4 - 

3.5 - 

3 - 

2.5 - 

2 - 

1.5 - 

1 - 

0.5 - 

- 



L=8, p, 

P2 

Pt 

P2 

i-1 



LP 




■HfiHBHmmiP 



0.5 



1.5 
T/J 



2.5 



FIG. 1 
tiplied 
T/J. 



The transverse and longitudinal spin stiffnesses mul- 
by L for L — 8 and L = 10 at p — 1 as a function of 



are built up by binning the values of the stiffnesses. 

Finite size scaling considerations for the stiffness 
p(T, L) predict that the singular part behaves as 



1 



p(T,L) = -f(L^t), 



(4) 



where t is the reduced temperature and v is the corre- 
lation length exponent. This scaling form predicts that 
plots of Lp(T, L) versus T should intersect at T c and that 
the L dependence of the slopes of the curves at this cross- 
ing point can be used to determine v. Fig. 1 shows the 
three principal stiffnesses multiplied by L as a function 
of T for sizes L = 8, 10. The longitudinal stiffness Lp^ 
exhibits a small peak which saturates with larger sizes. 
Both transverse stiffnesses behave identically and exhibit 
size dependent increases at low T proportional to L. Us- 
ing the scaling form in (4) for Lp\ and Lpi does not lead 
to a crossing point for different L. However, if we in- 
terpret the behavior of the longitudinal component Lp^ 
as a regular term and use the same scaling form for the 
difference Ap between the transverse and longitudinal 
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FIG. 2: The difference between transverse and longitudi- 
nal stiffnesses times L as a function of T/J for sizes L = 
4, 6, 8, ...20. The insert shows a finite size scaling collapse 
of the stiffnesses onto a single curve for T c = 1.4429(1) and 
i/ = 0.715(5). 



components, we obtain the results shown in Fig 2. The 
results clearly show a crossing temperature which moves 
to higher temperature as L is increased. 

Assuming that LAp(L, T) has the same phenomeno- 
logical behavior as the Binder cumulantii, we can de- 
termine the crossing point for different sizes L, L' and 
extrapolate T C (L) to infinite size in terms of the quantity 
(lnL'/I/) -1 where L is the smaller of the two sizes. For 
L = 8, 10 and 12 we find a good linear behavior for T C (L). 
Using this method, the critical temperature is estimated 
to be T c = 1.4429(1) and this value is in good agree- 
ment with the results from high-resolution Monte Carlo 
studiesi2ii^. The derivative of LAp with respect to t 
should behave as LAp' sw L 1 / i/ /'(0) and this prediction 
is obeyed quite well. The value of the static exponent v 
is estimated to be v — 0.715(5) which is remarkably close 
to the recently improved theoretical valued 4 . Using these 
values of T c and r^, we obtain a nice collapse of the data 
of Fig. 2 into a single universal curve as shown in the 
inset. 

Disordered Model 

The exponents obtained above for the pure Heisenberg 
model yield a negative value of a = 2 — dv ~ —.15. The 
Harris criterion 15 predicts that the critical behavior of 
the paramagnetic-ferromagnetic transition in the disor- 
dered system should have the same values if the specific 
heat exponent of the pure system is negative. Using the 
method described above, we have studied lattice sizes 
ranging from L = 4to_L = 16by averaging the principal 
values of the stiffness tensor over 500 (200) disorder con- 
figurations for small (large) sizes in the range 0.5 < p < 1. 
For p — .7 the finite size scaling analysis described above 
yields the values T c = 1.217(4), v — .70(1) and again al- 
lows for a nice collapse of the data with the pure Heisen- 
berg exponents. This behavior is in agreement with the 
predictions of the Harris criterion. 

Fig. 3 shows a schematic phase diagram for the model 
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FIG. 4: LAp plotted as a function of T/J for sizes L ranging 
from 8 to 16 for p = 0.6. The insets show a finite size scaling 
collapse of the data, (a) the upper transition has T c = 1.06(2) 
and v — 0.75(5) (b) the lower transition has Tn — 0.49(2) and 
v = 0.95(5). 




FIG. 3: Schematic phase diagram of the site disordered 
Heisenberg system. 



which differs from that reported by Matsubara et. al. 7 . 
We find both ordered and mixed phases but they do not 
meet at a decoupled tetracritical pointi At values of p 
slightly larger than 0.5, two transitions are observed. The 
upper transition T c is from the paramagnetic phase to a 
ferromagnetic phase where the transverse stiffnesses be- 
come different from the longitudinal stiffness. At a lower 
temperature, the two transverse components also acquire 
a difference which signals an order-order transition to a 
mixed phase. As indicated in Fig. 3, the upper transi- 
tion terminates close to p = 0.5 and there is a finite range 
of p where there is again a single transition to a mixed 
phase. We have studied the spin-spin structure factor of 
this phase and it has two sharp peaks corresponding to 
the two order parameters in equation (2). These order 
parameters are orthogonal to one another and hence the 
phase is characterized by a global chirality rh s x m and 
corresponds to a noncollinear order parameter. 

The stiffness results for p = 0.6 are shown in Figs. 



TABLE I: Results for the critical temperatures T c , Tm and 
the exponents v and (3 for various values of p. 
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4(a), (6), where it is clearly seen that two transitions 
are detected from the crossing of the stiffnesses Api3 
at T c ~ 1.06 and Ap 12 at T N ~ 0.49. In Fig. 4 (a) 
Api3 is the difference between the transverse and lon- 
gitudinal stiffnesses whereas in Fig. 4 (b) Api2 is the 
difference between the two transverse stiffnesses. These 
values of the critical temperatures agree with those re- 
ported previously^ for this same model at the equivalent 
concentration p — 0.4 but our value of the critical ex- 
ponent v at the upper transition is quite different. The 
previous work reported a value of v = 0.54 which cor- 
responds to a divergent specific heat. We find no such 
divergence and our value of v = 0.75(5) is compatible 
with the exponent of the pure ferromagnet as expected 
from the Harris criterion. Our value of v = 0.95(5) at the 
lower transition is somewhat larger than that reported by 
the previous authors. Hence the transition to the mixed 
" chiral" phase is characterized by a value of v ~ 1 which 
is close to the value reported for Heisenberg spin glasses 
and may belong to the same universality class. 

The values of the critical parameters obtained in the 
range 0.5 < p < 1 using the stiffness are given in Table I. 
The critical exponent (i is obtained from a finite size scal- 
ing of the order parameters in equation (2) and its value 
is fairly independent of p and is close to the pure Heisen- 



berg value. The transition to the mixed " chiral" phase at 
p = 0.5 where there is only a single transition is also de- 
scribed by a value of v ~ 1. Using v and /?, we estimate 
the value of the anomalous dimension exponent r\ ~ —0.3 
for all transitions to the mixed "chiral" phase. These val- 
ues of v and r\ are similar to those found in many spin 
glass mode ra 5 i 16 and may indicate that they are univer- 
sal. Although the model studied here has quenched site 
disorder, it may provide a useful reference model for the 
random bond models. By introducing dilution as well as 
competing interactions, we do not expect the correlation 
length exponent to change immediately since the mixed 
phase has a large negative value of a ~ — 1. However, for 
strong dilution, the order parameter could be much more 
complicated and the Harris argument may not apply. 

Our use of the helicity modulus in eqn(3) as a measure 
of the stiffness of the system assumes that the derivatives 
of the free energy are evaluated at zero twist which cor- 
responds to the equilibrium state in this model. For spin 
glasses^ and vortex glasses 1 ^, the equilibrium state does 
not necessarily correspond to zero twist and the helicity 
modulus would average to zero over the distribution of 
random bonds. The fluctuations of the helicity moduli 
would have to be studied in this case. 

In summary we have studied both the pure and site 
disordered 3d Heisenberg model on a simple cubic geome- 
try using Monte Carlo methods. We have shown that the 
spin stiffness provides a convenient method to determine 
transition temperatures as well as the correlation length 
exponent v. A mixed chiral phase which exhibits a uni- 
versal value of v ~ 1 is found at intermediate concentra- 
tions. Our results differ from those reported previously 
for the same model. 
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